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Abstract 

We study vortex clustering in type II Superconductors. We demonstrate that the "second 
peak" observed in magnetisation loops may be a dynamical effect associated with a density 
driven instability of the vortex system. At the microscopic level the instability shows up as 
the clustering of individual vortices at (rare) preferential regions of the pinning potential. 
In the limit of quasi-static ramping the instability is related to a phase transition in the 
equilibrium vortex system. 
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When the external magnetic field penetrates a type II superconductor, magnetic 
vortex lines appear inside the bulk of the sample. These vortex lines are repulsive and 
quantised each carrying a multiple q of the magnetic fiux quantum 0o. The energy per unit 
length of a vortex line is proportional to and accordingly one typically expects to have 
configurations with well separated single quantised vortex lines . We show below, however, 
that dynamical vortex clustering may become very important: vortices form inhomogeneous 
spatial structures which become relevant for the evolution of the system. This clustering 
may for instance cause the "second peak" observed in magnetisation loops and strongly 
affects the structure of the vortex system, magnetic relaxation, and the distribution of the 
local magnetic induction as measured in /i spin-relaxation experiments. The understanding 
of structural properties of vortices, underlying the presence of the "second peak" , is one of 
the central issues of current research in superconductivity, and is related to fundamental 
aspects of vortex matter ranging from dynamical behaviour to phase transitions, (see eg. 

Bi). 

In a superconducting sample in the presence of an external magnetic field, vortex 
lines penetrating from the surface into the bulk may be trapped on pinning centres Q 
leading to a spatially inhomogeneous vortex distribution and to a net magnetisation of the 
sample. As the external field is increased the vortex lines are squeezed together and therefore 
interact more strongly. Such a strong interaction will counteract the pinning forces and thus 
one expects, as in fact is usually observed, the magnetisation to decrease with increasing 
external field once vortices have fully penetrated the sample. 

However, upon further increase of the external field the magnetisation is often ob- 
served to increase again This behaviour leads to the so-called fishtail structure in the 
magnetisation data for YBCO or the arrowhead structure in equivalent data for BISCCO 
samples[^, 0. A similar peak structure has also been observed in low temperature super- 
conductors [Q. The second peak is one of the most important unsolved problems in vortex 
physics. It has been seen as a signature of a phase transition in the vortex system and its 
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explanation has often been attempted by focusing on collective aspects of pinning ||^, |^ or 
relaxation effects related to the sweeping of the field . Here we present evidence that 
dynamical effects are of crucial importance for the second peak, nevertheless the mechanism 
inducing the peak is, in the quasi-static limit, related to a density driven phase transition 
in the vortex system. 

In Fig. 1 we show for reference an example of half magnetisation loops in YBCO for 
a set of different temperatures. One sees how the maximum of the second peak moves to 
higher magnetic fields as the temperature is decreased. This suggests that the mechanism 
behind the increase in the magnetisation occurring after the first peak must also be active at 
zero-temperature. We demonstrate below that the origin of the second peak may be related 
to the possible grouping of vortices at favourable regions in the random pinning potential, 
leading to large local fluctuations in the vortex density. This "clustering" can occur in 
superconductors for which the ratio, n = A/^, between the magnetic penetration depth A 
and the coherence length ^ is not too small (from our Molecular Dynamics simulations we 
expect roughly n > 10). The clustering is a dynamical effect induced by the ramping of the 
external field, as is done in magnetisation experiments. As vortices entering (or leaving) 
the sample approach another vortex trapped at a position in the pinning potential three 
possibilities can occur. Firstly, the trapped vortex may be pushed ahead if the vortex- 
vortex repulsion is strong enough. Secondly, the approaching vortex may move around the 
trapped vortex. Or, if the trapped vortex is pinned by a force stronger than the maximum 
vortex-vortex repulsion, the approaching vortex may move into the favourable position in 
the pinning landscape in the immediate vicinity of the already trapped vortex. In this way 
the clustering may significantly enhance the effect of the rare strong pinning regions in an 
otherwise weak pinning background. 

For simplicity we here consider straight parallel flux lines at separation r for which 
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the interaction energy per unit length is 



U..{r) = ^ IMr/X') - Ko{r/0] , (1) 



where A' = A/ y/1 — b and ^' = ^/y2(l — b) are the effective field dependent London pen- 
etration depth and coherence length, respectively, b = B/Bc^ is the reduced magnetic 
induction relative to the upper critical induction i?c2, and Kq is a modified Bessel function. 
Vortices can be brought close together because of the attractive second term representing 
the interaction between the vortex cores. The maximum repulsive force calculated from 
Eq. (|l]) is an upper limit for the repulsive force between two vortex lines. The relative tilting 
and wiggling of vortex lines will lead to a significant decrease in the repulsion (see Eg. 0). 
Therefore, we expect the clustering effect to be at least as significant in three dimensions as 
we demonstrate here the effect to be in two dimensions. In fact clustering of vortices has 
been directly observed in electron microscopic imaging by Tonomura |Tn . 



In Fig. 2 we show the magnetisation obtained from zero temperature Molecular Dy- 
namics (MD) simulations. The vortex interaction is given by Eq.(P and it is cut off at half 
the system size. We use over damped dynamics in a square two dimensional system, with 
periodic boundary conditions, of sides 100.^. The external field is ramped by introducing 
vortices into a central strip with no pins the magnetisation plot is calculated by con- 
sidering the average density of vortices and its gradient in the pinned region. Fig. 2 clearly 
shows that an upturn in the magnetisation occurs as the external field is increased above 
the penetration field. The increase in the magnetisation coincides with the appearance of 
vortices clustered within areas of order It is important to mention that if the second 
term in Eq. |l| is left out[T^, (see dashed line in Fig. 2.), clustering of vortices cannot 



occur and no significant upturn in the magnetisation is observed in the simulations. At the 
highest field densities the potential in Eq.(|l]) loses its validity. This regime is also difficult 
to handle numerically. Hence we do not study the full half loop in the MD simulations but 
use below a simplified lattice model to study the increasing as well as the decreasing leg of 
the magnetisation. 



Let us now describe how the clustering of vortices prevents the magnetisation, M, 
from decaying with increasing field. We recall that approximately M oc jc, where jc is 
the critical current density produced by the volume pinning force Fp = Bjc- When a 
pinning centre becomes occupied by one or more trapped vortices the attractive short range 
pinning centre is effectively transformed into a longer range repulsive centre. Other vortices 
approaching this pin will feel a repulsive vortex- vortex force proportional to the number of 
trapped vortices at a distance A rather than the attractive force of range ^ from the initially 
"empty" pinning centre. The local pinning strength will fluctuate through the sample with 
typically a high density of weak pinning centres and only a few sparse local strong regions. 
Vortices clustered at the few strong pinning regions can then form spatially extended energy 
barriers which cage other diffusing vortices. We emphasise that this picture is somewhat 
schematic and that in reality dynamical and collective effects are important. Vortices are 
moving in an ever changing energy landscape produced by the combined effect of the static 
spatial pinning potential and the instantaneous metastable configuration of the interacting 
vortices. The importance of dynamical effects follows from the fact that the detailed form, 
especially the width, of the magnetisation loops dependence on the ramping rate of the 
external field. This is the case in experiments (see ref.s in ^ [13|) as well as in our 
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To study in detail the change of the effective energetic panorama due to clustering and 
its consequences for the "second peak" we consider now a schematic model. This approach 
is close in spirit to similar lattice systems introduced to describe fiuxons in superconductors, 
see e.g. In particular we study an extension of a coarse grained cellular-automaton-like 



model recently introduced by Bassler and Paczuski (BP) ||T6l. We consider a simplified 
version of a many body system with pair interactions given in eq. (|I|) representing a lattice 
model of repulsive particles in a pinning potential and in contact with a particle reservoir 
at a given density. Since our model explicitly allows multiple occupancy of lattice sites up 
to a value Nc2, we call it a Restricted Occupancy Model (ROM). We apply Monte Carlo 
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dynamics and use a Hamiltonian of the form: 

ij i i 

Here rij G {0, ...,Nc2} is an integer occupancy variable equal to the number of particles on 
site i. The parameter Nc2 plays the role of Bc2, it bounds the particle density per site below 
a critical value (here Nc2 = 27). The first term in eq.(0) represents the repulsive vortex 
interaction energy. The second term in eq.(|^) just normalises the particle self-interaction 
energy. Here we for simplicity choose the coarse graining length to be of order the zero 
temperature London penetration length A(0). This allows us to relate the restriction number 
Nc2 to the upper critical field Bc2 in the following way Nc2 = -Bc2A(O)^/0o where (po = hc/2e 
is the magnetic flux quantum. With this choice of coarse graining length it is natural as a 
first approximation to j are nearest neighbours 

and Aij = for all other couples of sites. We will below briefly discuss the validity of 
this approximation of the Aij matrix. The third term in eq.(|^) represents a random pinning 
potential acting on a fraction p = 0.5 of the lattice sites with = 0.5 and A^ = elsewhere. 
(The same set of interactions is used in the BP model [16|). Two opposite sides of our square 



system (L = 32 ) are in contact with a reservoir at a given density Ng^t- Particles are 



introduced and escape the system through the reservoir only. 

Fig. 3 shows the results of our Monte Carlo simulations of this model. We ramp 
Next and record the magnetisation, M = — Next (with = ^i) / L"^) , as a function 
of Next- The ramping of the external reservoir density, Next, is simply done by increasing 
it from zero up to some given value (and then decreasing back) by a sequence of small 
increments, ANq. After each increment the system is let to relax for a time r (in unit of 
Monte Carlo sweeps). This corresponds to a sweep rate of the apphed field of 7 = ANq/t. 

We recorded magnetisation loops for several values of the ratio k* = (\nAy/An)~^ 
which controls the interaction potential. The parameter k* is qualitatively similar to k = 
X/^. When k* is large enough a definite second peak appears in M. These magnetisation 
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loops looks qualitatively similar to the experimental loops, though the peak position and 
amplitude of peaks in Fig. 3 are much more asymmetric than the peaks of the experimental 
data in Fig. 1. Nevertheless, typically magnetisation loops in other experimental samples 
and compounds are clearly asymmetric, see e.g. In fact even Fig. 1 exhibits, on 

careful inspection, slight asymmetries. 

The detailed loop shape depends not only on the choice of k*, as shown in Fig. 3, 
but also on the value of Nc2 and the ramping rate, 7. We find that in the 7 — > limit, the 
second peak location, is associated with a sharp jump in Mgg = lim^^o ^(7); corresponding 



to a true transition ||Tj]. The precise nature of this transition is currently under study. The 
transition occurs above the melting transition and leads to a significant increase in the 
effective energy barriers experienced by the diffusing vortices. 

Increasing Nc2 (corresponding to higher values of Bc2) increases the separation be- 
tween the first and the second peak. The specific features of the loops do also depend on the 
ratio between the characteristic relaxation time of the vortex system and the ramping rate. 



for details see ||14| , [181 . Quantitative differences between the simulations and experiments 
are to be expected. One reason is that the interaction strength between vortex lines depends 
on the magnetic induction and the temperature. As for instance in the London approxima- 
tion of Eq.(|l]) through the field and temperature dependence of A' and ^' . This corresponds 
to a field and temperature dependence of the individual elements of the coupling matrix 
Aij in Eq.(H). Moreover, the effective vortex screening length. A' (see Eq. |l]), increases 
with increasing magnetic field. This effect implies that non-zero elements between sites 
of separation larger than nearest neighbour may become relevant as the field is increased. 
Interestingly, however, the present simple approximation captures the qualitative features 
of the magnetic properties. 

We emphasise that according to the above picture the second peak is a dynamical 
effect associated with the vortices being forced in and out of the sample. It is linked 
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to a density driven instability of the vortex system. At the more microscopic level of 
description used in the Molecular Dynamics simulations the instability shows up as the 
clustering of individual vortices at (rare) preferential regions of the pinning potential. In 
the coarse grained description of the considered lattice model the instability is related to 
an underlying density driven phase transition of the equilibrium system. The instability 
induces a dramatic change of the effective collective energy landscape encountered by the 
diffusing vortices. This, in turn, enormously enhances equilibrium times and induces 



the presence of significant more spatial disorder, in strict correspondence with glass formers 
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This clustering or density instability can produce the second peak in the magneti- 
sation loops. It can be thought of as a type of strong plastic deformation, an effect which 
should be observable in neutron scattering or /i-spin resonance experiments probing the 
distribution of local magnetic induction. 
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Figures Captions 



Fig. 1 An example of magnetisation data and the fishtail effect in a single crystal of 
YBa2Cu408 at temperautres 20K 50K and 65K. The magnetic moment is measured con- 
tinuously using a vibrating sample magnetometer while the applied field is ramped from 
zero to eight Tesla and back to zero at a rate lOmT/sec . The 124 system is neither sus- 
ceptible to twinning or oxygen inhomogeneity (unlike YBa2Cu307-d) indicating that the 
fishtail effect is intrinsic to the pinning of point-like disorder in the crystal lattice. 

Fig. 2 Magnetisation (dB/dx) vs field (B) for simulation with k, of 67. The external field 
is ramped by adding one vortex between relaxation intervals of 40 where each time step 
is a maximum of 0.01. Fields calculated for ^ of 15 A. Solid line represents the system 
interacting through Eq. 1 (soft core). The dashed fine represents the system interacting 
through Eq. 1 with last term omitted (hard core). The short-long dashed curve represents 
the proportion of stacked vortices in the soft core case. 

Pinning centres are represented by Gaussian wells of width ^ and of amplitude 0.3^^ 
times the condensation energy. In this simulation they can exert a maximum pinning force 
(at zero external field) of 4 • IO^^A/ctv? for the k = 10 and 8 ■ lO^A/cim? for n = 100. 

Fig. 3 The magnetisation, M, as a function of the applied field density, A^^..^.^, in the 2D 
R.O.M. model for k* = 0.43, 0.76, 0.79. The ramp rate for Next is ANo/t = 10"^. 
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